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The effects of non-zero baryonic chemical potential on the drag and diffusion coefficients of heavy 
quarks propagating through a baryon rich quark gluon plasma have been studied. The nuclear 
suppression factor, Raa for non-photonic single electron spectra resulting from the semileptonic 
decays of hadrons containing heavy ffavours have been evaluated for low energy collisions. The role 
of non-zero baryonic chemical potential on -Raa has been highlighted. 
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I. INTRODUCTION 



The nuclear collisions at low energy RHIC run [l|, 
H and GSI-FAIR is expected to create a ther- 
mal medium with large baryonic chemical potential 
(Pb) and moderate temperature (T). The heavy 
flavours namely, charm and bottom quarks may play 
a crucial role in understanding the properties of such 
medium because they do not constitute the bulk 
part of the system and their thermalization time 
scale is larger than the light quarks and gluons and 
hence can retain the interaction history very effec- 
tively. The perturbative QCD (pQCD) calculations 
indicate that the heavy quark (Q) thermalization 
time, is larger 0, Q than the light quarks and 
gluons thermalization scale r^. Gluons may ther- 
malized even before up and down quarks 0| • In 
the present work we assume that the Quark Gluon 
Plasma (QGP) is formed at time r^. Therefore, 
the interaction of the non-equilibrated heavy quarks 
with the equilibrated QGP for the time interval 
Ti < T < can be treated within the ambit of the 
Fokker-Planck (FP) equation i.e. the heavy 

quark can be thought of executing Brownian mo- 
tion 0, H, [l3-[i3 in ttie heat bath of QGP during 
the said interval of time. 

As the relaxation time for heavy quarks of mass 
M at a temperature T are larger than the corre- 
sponding quantities for light partons by a factor of 
M/T{> 1) [3| i.e. the light quarks and the gluons 
get thermalized faster than the heavy quarks, the 
propagation of heavy quarks through QGP (mainly 
contains light quarks and gluons) therefore, may be 
treated as the interactions between equilibrium and 
non-equilibrium degrees of freedom. The FP equa- 
tion provide an appropriate frame work for such pro- 
cesses. In case of low energy collisions the radiative 
energy loss of of heavy quarks will be much smaller 
than the loss due to elastic processes. Moreover, the 
thermal production of charm and bottom quarks can 
be ignored for the range of temperature and bary- 
onic chemical potential under study. Therefore, the 



FP equation is better applicable in the present situ- 
ation. 

The paper is organized as follows. In the next 
section a brief description of Fokker Planck equa- 
tion and the T and fj,q dependence of the drag and 
diffusion co-efficients are outlined, the non-photonic 
electron spectra is discussed in section III, the initial 
conditions and the space time evolution have been 
discussed in section IV, section V is devoted to the 
discussions on nuclear suppression and finally sec- 
tion VI contains the summary and conclusions. 



II. THE FOKKER PLANCK EQUATION 

The Boltzmann transport equation describing a 
non-equilibrium statistical system reads: 
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where p and E denote momentum and energy, Vx 
(Vp) are spatial (momentum space) gradient and 
f{x,p,t) is the phase space distribution (in the 
present case / stands for heavy quark distribution). 
The assumption of uniformity in the plasma and ab- 
sence of any external force leads to 



dl 
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The collision term on the right hand side of the above 
equation can be approximated as (see [HI, for 
details) : 
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where we have defined the kernels 



^ J d?kuj{p, k)ki 
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D5ij where 



for I p 1^ 0, ^ 7Pi and B^j 
7 and D stand for drag and diffusion co-efficients 
respectively. The function a;(p, k) is given by 



i2ny 
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where /' is the phase space distribution, in the 
present case it stands for light quarks and gluons, 
V is the relative velocity between the two collision 
partners, a denotes the cross section and g is the 
statistical degeneracy. The co-efficients in the first 
two terms of the expansion in Eq. [3] are compara- 
ble in magnitude because the averaging of ki in- 
volves greater cancellation than the averaging of the 
quadratic term kikj. The higher power of fc^'s are 
smaller [8]. 

With these approximations the Boltzmann equa- 
tion reduces to a non-linear integro-differential equa- 
tion known as Landau kinetic equation: 
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The nonlinearity is caused due to the appearance of 
/' in Ai and Bij through w{p, k). It arises from the 
simple fact that we are studying a collision process 
which involves two particles - it should, therefore, 
depend on the states of the two participating parti- 
cles in the collision process and hence on the prod- 
uct of the two distribution functions. Considerable 
simplicity may be achieved by replacing the distri- 
bution functions of one of the collision partners by 
their equilibrium Fermi-Dirac or Bose-Einstein dis- 
tributions (depending on the statistical nature) in 
the expressions of Ai and Bij . Then Eq. |6] reduces 
to a linear partial differential equation - usually re- 
ferred to as the Fokker-Planck equation[9] describing 
the interaction of a particle which is out of thermal 
equilibrium with the particles in a thermal bath of 
light quarks, anti-quarks and gluons. The quantities 
Ai and Bij are related to the usual drag and diffu- 
sion coefficients and we denote them by 7^ and Dij 
respectively {i.e. these quantities can be obtained 
from the expressions for Ai and Bij by replacing 
the distribution functions by their thermal counter- 
parts): 

The evolution of the heavy quark distribution (/) 
is governed by the FP equation 11] 
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where ji and Dij are the drag and diffusion coef- 
ficients. The elastic collisions of heavy quarks {Q 
stands for charm and bottom) with thermal light 
quarks (q), anti-quarks (q) and gluons (g) i.e. Qq 
Qq, Qq — ^ Qq and Qg Qg have been used to eval- 
uate the drag and diffusion coefficients as indicated 



in [llj, [l8| . The thermal distribution of the quarks 
and the anti-quarks are responsible for the T and fj,q 
dependence of the drag and diffusion co-efficients. 
The gluon distribution introduces the T dependence 
to these quantities, chemical potential dependence 
of the drag and diffusion co-efficients originate from 
the thermal phase space of quarks and anti-quarks. 



A. The drag and diffusion co-efRcients 



At low y/sNN the net baryon density at mid- 
rapidity is non-zero and its value could be high de- 
pending on the value of y'SNAr- Therefore, we need 
to solve the FP equation for non-zero fj,B- The drag 
and diffusion coefficients are functions of both the 
thermodynamical variables: fj,B and T. 

The energy dependence of the chemical potential 
has been obtained from the parametrization of the 
experimental data on hadronic ratios as pJi] (see 
also 0), 



(8) 



where a = 0.967 ± 0.032 GeV and b = 6.138 ± 0.399 
GeV. The parametrization in Eq. |5] gives the values 
of /is at the freeze-out. The corresponding values 
at the initial condition are obtained from the baryon 
number conservation equation. The initial baryonic 
chemical potential carried by the quarks fJ-q{— Ms/3) 
are shown in table 1 for various ^sntv under consid- 
eration. 

In the present work we take as = 0.3 because the 
dependence of the strong coupling on the temper- 
ature and baryonic chemical potential is not accu- 
rately known yet. The sensitivity of collisional en- 
ergy loss on the running is studied in Ref. [2lj 
in detail. The variation of the drag coefficients of 
charm quarks (due to its interactions with quarks 
and anti-quarks) on the baryonic chemical potential 
for different T are displayed in Fig. [T] The drag co- 
efficient for the process : Qg — )■ Qg is ~ 8.42 x 10^'^ 
fm-i (1.86 X 10-2 fm-i) for T = 140 MeV (190 
MeV) (not displayed in Fig. [T]). The T and fig de- 
pendence of the drag and diffusion co-efRcients may 
be understood as follows. The drag may be defined 
as the thermal average of the square of the invari- 
ant transition amplitude weighted by the momentum 
transfer for the reactions qQ — > qQ, Qq — Qq and 
gQ — > gQ. As the temperature of the thermal bath 
increases the light quarks (q) and the gluons move 
faster and gain the ability to transfer larger momen- 
tum during their interaction with the heavy quarks 
- resulting in the increase of the drag of the heavy 
quarks propagating through the partonic medium. 
Since the average momentum of the quarks increases 
with fiq, similar behaviour is expected in the varia- 
tion of drag with baryonic chemical potential. This 
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trend is clearly observed in the results displayed in 
Fig. [T] for charm quark. The drag due to the pro- 
cess Qq — >■ Qq is larger than the Qq — >■ Qq interac- 
tion because for non-zero chemical potential, the Q 
propagating through the medium encounters more q 
than g at a given /i^. For vanishing chemical poten- 
tial the contributions from quarks and anti-quarks 
are same. 

In the same way it may be argued that the dif- 
fusion coefficient involves the square of the momen- 
tum transfer - which should also increase with T and 
jjLq as observed in Fig. [51 The diffusion co-efficient 
for charm quarks due to its interaction with gluons 
is given by - 1.42 x lO'^ GeVVfm (4.31 x lO^^ 
GeVVfm) for T 140 MeV (190 MeV). It may 
be mentioned here that the drag increases with T 
when the system behave like a gas. In case of liq- 
uid the drag may decrease with temperature (ex- 
cept very few cases) - because a substantial part of 
the thermal energy goes in making the attraction 
between the interacting particles weaker - allowing 
them to move more freely and hence making the drag 
force lesser. The drag co-efficient of the partonic 
system with non-perturbative effects may decrease 
with temperature as shown in Ref. 22| - because in 



this case the system interacts strongly more like a 
liquid. The heavy quark momentum diffusion co- 
efficient has been computed [l^ at next to leading 
order within the ambit of hard thermal loop approxi- 
mations. For T ^ 400 MeV the momentum averaged 
pQCD value (for jiq — 0) of the diffusion co-efficient 
obtained in the present work is comparable to the 
value obtained in [23|] in the leading order approxi- 
mation for the same set of inputs {e.g. strong cou- 
pling constant, number of flavours etc). The drag 
and diffusion coefficients for bottom quarks are dis- 
played in Figs. [3] and 2] respectively, showing qual- 
itatively similar behaviour as charm quarks. The 
drag co-efficients for bottom quarks due to the pro- 
cess Qg — Qg is given by ~ 3.15 x lO^'^ fm^^ and 
6.93 X 10-3 fm-i at T = 140 McV and 190 MeV re- 
spectively. The corresponding diffusion coefficients 
are - 1.79x10^3 GeVVfm and 5.38x10-^ GeVVfm 
at T=140 MeV and 190 MeV respectively. 



gg QQ and qq QQ. Here we intend to deal 
with the nuclear suppression factor, Raa, which in- 
volves the ratio of the momentum distribution func- 
tions. Therefore, the final results may not be too 
sensitive to the initial distributions because of some 
cancellations that may take place in the ratio. 
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FIG. 1: Variation of the drag coefficient of charm quark 
due to its interactions with light quarks and anti-quarks 
as a function of /iq for different temperatures. 
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III. THE NON-PHOTONIC ELECTRON 
SPECTRA 

After obtaining the drag and diffusion coefficients 
we need the initial heavy quark momentum distri- 
butions for solving the FP equation. For low col- 
lision energy rigorous QCD based calculations for 
heavy flavour production is not available (for higher 
y/sWN = 200 GeV see Ref. [H for rigorous QCD 
calculations). In the present work this is obtained 
from pQCD calculation [l^, for the processes: 



FIG. 2: Variation of the diffusion coefficient of charm 
quark due to its interactions with light quarks and anti- 
quarks as a function of /i, for different temperatures. 



With the initial condition mentioned above the 
FP equation has been solved for the heavy quarks. 
We convolute the solution with the fragmenta- 
tion functions of the heavy quarks to obtain the 
Pt distribution of the heavy (B and D) mesons 
{dN^^^/qTdqr). For heavy quark fragmentation we 
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use Peterson function [27| given by: 

1 



(9) 



for charm quark ec = 0.05. For bottom quark Cb = 
{Mc/Mbfec where [Mb] is the charm (bottom) 
quark mass. 
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FIG. 3: Same as Fig. [T] for bottom quark. 
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FIG. 4: Same as Fig. [2] for bottom quark. 



The non-photonic single electron spectra originate 
from the decays of heavy flavour mesons - e.g. D 
Xev at mid-rapidity (y = 0) can be obtained as 

follows [ii-Ii^: 

= I dqT—^F{pT,qT) (10) 



qrdqj 



where 



F{pT. qT)^oj f ^(Pt-Qt) g(p^ q^/M) (11) 
J 2pTPTqT 



where M is the mass of the heavy mesons (D or B), 
uj = 96(1 - 8m2 + Sm^ - - l2mHnw?)-^M-^ 
(m — Mx/M) and g{Ee) is given by 



9{Ee 



(M - 2E,) 



(12) 



related to the rest frame spectrum for the decay D 
Xev through the following relation [28] 



1 dV 



ujg{Ee). 
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We evaluate the electron spectra from the decays 
of heavy mesons originating from the fragmentation 
of the heavy quarks propagating through the QGP 
medium formed in heavy ion collisions. In the same 
way the electron spectrum from the p-p collisions 
can be obtained from the charm and bottom quark 
distribution which goes as initial conditions to the 
solution of FP equation. The ratio of these two 
quantities gives the nuclear suppression, Raa as : 
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called the nuclear suppression factor, will be unity in 
the absence of any medium. In the above equation 
A'coU denotes the number of nucleon nucleon colli- 
sions in Au-|-Au interaction. However, the experi- 
mental data [m, im at RHIC energy {y/sNN=200 
GeV) shows substantial suppression {Raa < 1) for 
> 2 GeV indicating substantial interaction of 
the plasma particles with charm and bottom quarks 
from which electrons are originated through the 
process: c{b) (hadronization) — >■ D(i?) (decay) — > 
e + X. The loss of energy of high momentum heavy 
quarks propagating through the medium created in 
Au-I-Au collisions causes a depletion of high px elec- 
trons. 



IV. THE INITIAL CONDITIONS FOR THE 
SPACE-TIME EVOLUTION 

The nuclear suppression for heavy quarks depend 
on the parameters like initial temperature (T^), ther- 
malization time (r^), equation of state (EOS) and 
the transition temperature (Tc). 

We assume that the system reaches equilibration 
at a time Ti after the collision at temperature Ti 
which are related to the produced hadronic (predom- 
inantly mesons) multiplicity through the following 
relation: 
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FIG. 5: The nuclear suppression factor Raa as a func- 
tion of pt due to the interaction of the charm quark 
(solid line) and anti-quark (dashed-dot line) for /i, — 200 
MeV. The net suppressions including the interaction 
of quarks, anti-quarks and gluons for /i, = 200 MeV 
(dashed hne) and Hq — (with asterisk) are also shown. 



where Ra is the radius of the system, ({3) is the 
Riemann zeta function and Oe// — ^^5e///90 where 
g^fj (= 2 X 8 -h 7 X 2 X 2 X 3 X is the degen- 

eracy of quarks and gluons in QGP, NF='DximheT of 
flavours. 

The value of the muhiplicities for various -y/sR/v 
have been calculated from the Eq. below [lH ; 
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Ncoii is the number of collisions and contribute x 
fraction to the multiplicity dupp/dy measured in pp 
collision. The number of participants, Npart con- 
tributes a fraction (1 — x) to dripp/dy, which is given 

by 



dn^ 



VP 



dy 



^2.5-0.25ln{s) + 0m3ln^{s) (17) 



The values of Npart and Ncoii are estimated for 
(0 — 5%) centralities by using Glauber Model. The 
value of X depends very weakly on ^/stm [13; in 
the present work we have taken x = 0.1 for all the 
energies. 

The time evolution of the temperature and the 
baryon density (ub) have been obtained by solving 
the following equations: 



0, 







(18) 



in (1+1) dimension with boost invariance along the 
longitudinal direction [3^1 • In the above equation 



TABLE I: The values center of mass energy , dN/dy, 
initial temperature (Ti) and quark chemical potential - 



used in the present calculations. 



V(sNiv)(GeV) 


dJV 
dy 


r,(MeV) /i,(MeV) 


39 


617 


240 62 


27 


592 


199 70 


17.3 


574 


198 100 


7.7 


561 


197 165 



T'"^ = {e + P)u^u'^ — g^'^P, is the energy momen- 
tum tensor and rig — ubu^ is the baryonic flux, 
where e is the energy density, P is the pressure, 
and is the hydrodynamic four velocity. The ra- 
dial co-ordinate dependence of T and n b have been 
parametrized as in Ref. [1^. The velocity of sound 
for the QGP phase is taken as Cs = l/-\/4. Some 
comments on the effects of the radial flow are in 
order here. The radial expansion will increase the 
size of the system and hence decrease the density of 
the medium. Therefore, with radial flow the heavy 
quark will traverse a larger path length in a medium 
of reduced density. These two oppositely competing 
effects may have negligible effects on the nuclear sup- 
pression (see also [i3|)- Moreover, at lower collision 
energies (as in the present case) the amount of ra- 
dial flow will not be as substantial as in -^snIv = 200 
GeV. 

The total amount of energy dissipated by a heavy 
quark in the QGP depends on the path length it tra- 
verses. Each parton traverses different path length 
which depends on the geometry of the system and 
on the point where its is produced. The probabil- 
ity that a parton is created at a point (r, 0) in the 
plasma depends on the number of binary collisions 
at that point which can be taken as [l3l |: 



(19) 



where R is the nuclear radius. A parton created at 
(r, (p) in the transverse plane propagate a distance 
L — R^ — r^sin^cj) — rcoscf) in the medium. In the 
present work we use the following equation for the 
averaging of the drag coefficient : 



rdrd(j)P{r, < 



L/v 



dTj{T) (20) 



where v is the velocity of the propagating partons. 
Similar averaging has been performed for the dif- 
fusion co-efficient. For a static system the T and 
liq dependence of the drag and diffusion co-efficients 
of the heavy quarks enter via the thermal distribu- 
tions of light quarks, anti-quarks and gluons through 
which it is propagating. However, in the present 
scenario the variation of the temperature and the 
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baryon density with time are governed by the equa- 
tion of state or the velocity of sound of the thermal- 
ized system undergoing hydrodynamic expansion. In 
such a scenario the quantities hke F (Eq. [20)1 and 
hence Raa becomes sensitive to velocity of sound in 
the medium. 
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FIG. 6: Nuclear suppression factor, -Raa as function of 
PT for various ^snn- Inset: the nuclear suppression 
factor due to the interaction of D meson in a thermal 
medium of pions and nucleons. 



V. THE NUCLEAR SUPPRESSION 

To demonstrate the effect of non-zero baryonic 
chemical potential we evaluate Raa for fiq = 200 
MeV and = for a given T, = 200 MeV. The re- 
sults are displayed in Fig. [5] - representing the com- 
bined effects of temperature and baryon density on 
the viscous drag and diffusion. The viscous drag on 
the heavy quarks due to its interaction with quarks 
is larger than that of its interactions with the anti- 
quarks (FiglD). Resulting in larger suppression in 
the former case than the later. The net suppression 
of the electron spectra from the Au-|-Au collisions 
compared to p-l-p collisions is effected by quarks, 
anti-quarks and gluons. The results for net sup- 
pressions are displayed for fiq ~ 200 MeV (dashed 
line) and Hq = (with asterisk). The experimental 
detection of the non-zero baryonic effects will shed 
light on the net baryon density (and hence baryon 
stopping) in the central rapidity region. However, 
whether the effects of non-zero baryonic chemical po- 
tential is detectable or not will depend on the overall 
experimental performance. 

The results for Raa are shown in Fig. [6] for vari- 
ous y/SNN with inputs from table I. We observe that 



at large pt the suppression is similar for all energies 
under consideration. This is because the collisions 
at high are associated with large temperature 

but small baryon density at mid-rapidity- which is 
compensated by large baryon density and small tem- 
perature at low ^ysNN collisions. Low pT particles 
predominantly originate from low temperature and 
low density part of the evolution where drag is less 
and so is the nuclear suppression. 

So far we have discussed the suppression of the 
non-photonic electron produced in nuclear collisions 
due to the propagation of the heavy quark in the the 
partonic medium in the pre-hadronization era. How- 
ever, the suppression of the D mesons in the post 
hadronization era (when both the temperature and 
density are lower than the partonic phase) should in 
principle be also taken into account. We have esti- 
mated Ra a for D mesons due to ts interaction with 
pions |36| and nucleons and found that it has a 
value value closer to unity, indicating the fact that 
the hadronic medium (of pions and nucleons) is un- 
able to drag the D mesons. Therefore, the measured 
depletion in Raa for the non-photonic electron will 
indicate the presence of partonic medium and the 
amount of depletion may the used to characterize 
the thermal medium. 

It has been shown in [s^ that a large enhance- 
ment of the pQCD cross section is required for the 
reproduction of experimental data on elliptic flow at 
RHIC energies. In our earlier work fis'l we have eval- 
uated the Raa for non-photonic single electron spec- 
tra resulting from the semileptonic decays of hadrons 
containing heavy flavours and observed that the data 
from RHIC collisions at ^/SNN = 200 GeV are well 
reproduced by enhancing the pQCD cross sections 
by a factor 2 and with an equation of state P = e/4. 
In the same spirit we evaluate Raa with twice en- 
hanced pQCD cross section and keeping all other 
quantities unaltered (Fig. [7]). The results in Fig. [7] 
show stronger suppression as compared to the re- 
sults displayed in Fig. [6l but it is similar in all the 
energies under consideration. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the effects of baryonic chemical 
potential and temperature on the drag and diffusion 
coefficients of heavy quarks moving in a thermalized 
system of quarks and gluons. We have observed that 
both the drag and diffusion coefficients increase with 
temperature and chemical potential. When we have 
enhanced the pQCD cross section for the interac- 
tion of the heavy quarks with the thermal system 
by a factor of two - the resulting suppressions in 
Raa are between 20% - 30% iorJJ^ = 39 - 7.7 



GeV. The radiative energy loss |39l443| (see [4J| for 
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FIG. 7: Same as Fig.|6]with enhancement of cross section 
by a factor of 2. 



a review) of heavy quarks is suppressed due to dead 
cone effects and has been neglected in the present 



work. Moreover, at low collisions energies the coUi- 
sional loss [2l|, |4^, 46 1 is dominant over its radiative 
counter part (see 47| for details). It may be men- 



tioned here that the theoretical formalism, the FP 
equation is applicable better for heavy quarks than 
light quarks and gluons (because of their frequent 
productions and annihilations). However, the pro- 
duction of charm and bottom quarks are smaller at 
low energy collisions making the measurements of 
non-photonic single electron spectra and hence Raa 
for heavy quarks difficult. The detection of the non- 
zero baryonic chemical potential effects observed in 
the present work through the nuclear suppression 
factor will help in determining the net baryon den- 
sity (and hence baryon stopping) in the mid-rapidity 
region. However, whether such effects is detectable 
experimentally or not will depend on the overall ex- 
perimental performance. 
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